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Abstract 

A  model  for  a  100  kW  class  solid  oxide  fuel  cell  (SOFC)  system  running  on  biogas  from  a  sewage  sludge  digestion  plant  was  implemented 
in  a  process  flow  scheme  using  external  steam  reforming.  The  model  stack  consisted  of  planar  anode  supported  cells  operated  at  800  °C 
displaying  state-of-the-art  electrochemical  performance  (0.15  W/cm2  at  80%  fuel  utilisation).  Real  annual  data  from  an  existing  sewage 
plant  were  used  as  input  to  the  model.  From  the  input  of  43  m3/h  biogas  (63%  CH4),  equivalent  to  269  kW  (higher  heating  value,  HHV), 
the  SOFC  stack  was  calculated  to  deliver  131  kWei  electricity  (48.7%)  using  a  steam-to-carbon  ratio  of  0.5.  This  would  allow  the  sewage 
site  to  more  than  cover  its  own  electrical  needs,  hence  to  depollute  the  waste  stream  at  negative  energy  cost.  In  its  current  exploitation 
using  a  low  efficient  gas  engine  (130  kW),  the  site  is  only  ~50%  self-sufficient.  Special  attention  was  given  to  the  thermal  balance  of  the 
stack.  The  stack  developed  heat  (143  kW)  could  be  balanced  by  endothermal  reforming  (78  kW)  and  by  cathode  excess  air  X  (=3),  allowing 
a  temperature  difference  between  stack  inlet  and  outlet  of  200  K.  The  case  was  compared  to  other  fuel  scenarios.  Steam-added  biogas 
behaves  basically  identically  to  steam-reformed  methane.  For  partial  oxidation  of  biogas  or  pure  hydrogen  feeding,  electrical  efficiency 
drops  to  under  43%  while  X  needs  to  be  raised  to  4.5  to  maintain  the  200  K  thermal  gradient  over  the  stack. 
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1.  Introduction 

Conversion  of  biogas  to  electricity  presents  an  attractive 
niche  application  for  fuel  cells  [1].  Biogas  production  sites 
are  small  (of  a  few  kWei  to  few  MWei),  plentiful  (from 
sewage  sludge,  farm  waste,  landfill,  industrial  liquid  waste, 
municipal  and  industrial  organic  solids)  and  steady  in  sup¬ 
ply.  Biogas  produced  from  sewage  sludge  digestion  is  par¬ 
ticularly  useful  because  of  its  stable  composition,  containing 
a  high  fraction  of  methane  and  a  low  level  of  sulphur. 

Electricity  production  on  average  sewage  treatment  plant 
sites  (about  100kWei)  is  typically  obtained  from  gas  en¬ 
gines.  However,  depending  on  lifetime  and  operating  con¬ 
ditions  of  the  engine,  electrical  conversion  is  then  relatively 
low  [1]  and  may  not  suffice  to  cover  the  site’s  annual  elec¬ 
tricity  needs,  while  heat  production  recovered  from  the  en¬ 
gine  can  be  in  excess  of  the  site  requirements,  i.e.  space 
heating  and  maintaining  the  digesters  temperature  (35  °C). 
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In  such  a  situation,  the  electricity  deficit  has  to  be  bought 
in  and  the  excess  heat  rejected.  Fuel  cells,  owing  to  higher 
electrical  efficiency,  could  thus  better  fit  the  plant  needs  with 
an  output  higher  in  electricity  and  lower  in  heat,  from  the 
on-site  generated  biogas. 

Amongst  fuel  cells  for  biogas  conversion,  the  high  temper¬ 
ature  solid  oxide  fuel  cell  (SOFC)  type  is  especially  suited 
because  of  the  capability  of  thermally  integrated  biogas  re¬ 
forming  and  because  of  manageable  tolerance  against  fuel 
contaminants  [1].  A  1  kWei  SOFC  unit  (900  °C)  from  Sulzer 
HEXIS,  Switzerland,  was  successfully  operated  (28  %  lower 
heating  value  (FHV)  efficiency)  on  farm  biogas  in  Switzer¬ 
land  for  1  year  [2].  Fuel  Cell  Technologies  (FCT)  from 
Canada  demonstrate  5kWei  SOFC  units  (46%  FHV  effi¬ 
ciency),  also  for  biogas  application  [3].  A  tubular  SOFC 
system  of  100kWei  (1000  °C)  has  been  tested,  not  with  bio¬ 
gas,  but  with  pipeline  natural  gas,  for  a  successful  2-year 
run  [4]  achieving  47%  FHV  electrical  efficiency.  Because 
of  the  similarities  between  natural  gas  and  sewage  biogas 
in  fuel  reforming  and  final  sulphur  levels,  such  a  system 
can  no  doubt  be  operated  also  on  the  latter  of  the  two 
fuels.  The  same  argument  applies  to  a  planar  SOFC  sys- 


128 


J.  Van  herle  et  al.  / Journal  of  Power  Sources  131  (2004)  127-141 


tern,  shown  recently  to  achieve  5.4kWei  (38%  efficiency) 
on  steam-reformed  methane  [5].  This  unit  operated  at  lower 
temperature  (800  °C)  owing  to  the  use  of  anode  supported 
thin  electrolyte  cell  (ASE)  fabrication  technology  (20  cm  x 
20  cm  cells,  361  cm2  active  area),  that  possess  lower  inter¬ 
nal  resistance  loss.  This  technology  (planar  ASE)  will  be 
referred  to  in  the  system  model  study  presented  here. 

A  previous  system  model  study,  for  an  SOFC  unit  of  a 
few  kWei  converting  farm  biogas  using  partial  oxidation  as 
fuel  processing  step,  was  published  by  the  present  authors 
[6] .  Partial  oxidation  had  been  chosen  as  the  most  practical 
and  economical  solution  for  such  a  small  system.  Phenom¬ 
ena,  such  as  gas  diffusion  limitation  and  the  risk  of  carbon 
deposition,  left  unquantified  in  the  previously  used  model, 
were  implemented  in  the  presently  improved  model. 

This  study  is  concerned  with  a  sewage  treatment  plant, 
producing  roughly  100  kWei -equivalent.  At  this  size,  steam 
reforming  as  fuel  processing  step  is  preferred  to  partial  oxi¬ 
dation,  for  its  higher  efficiency  and  more  appropriate  thermal 
integration  (endothermicity)  with  the  fuel  cell  stack.  Special 
attention  to  thermal  management  will  be  given  here.  The 
different  fuel  processing  options  will  be  compared. 

Real  data  from  an  existing  sewage  sludge  treatment  plant, 
producing  continuously  approximately  1000  m3  biogas/day, 
are  taken  as  input  to  the  SOFC  system  model.  The  aim  is  to 
demonstrate,  using  the  SOFC  instead  of  the  gas  engine,  that 
site  self-sufficiency  in  electricity  and  heat  production  from 
the  generated  biogas  can  be  attained,  resulting  in  depollution 
of  its  waste  stream — which  constitutes  the  primary  purpose 
of  the  plant — at  no  net  energy  cost. 

2.  Experimental 

Biogas  samples  from  sewage  sludge  digestion  plants  (two 
different  sites)  have  been  periodically  collected  in  glass  bot¬ 
tles  and  analysed  for  their  chemical  composition,  using  a  se¬ 
ries  of  techniques.  The  levels  of  the  main  constituents  were 
determined  by  gas  chromatography  (GC),  using  a  Micro-Gas 
Chromatograph  CP2300  from  Varian  Inc.,  Zug,  Switzerland, 
equipped  with  a  Poraplot-Q  column  and  thermal  conductiv¬ 
ity  detector.  The  averaged  composition  was  used  as  input  to 
the  SOFC  system  model,  together  with  site  production  data, 
such  as  the  hourly  biogas  flow  rate. 

Microcontaminant  levels  were  determined  by  total  com¬ 
bustion  with  O2  followed  by  ion-exchange  chromatogra¬ 
phy  (for  Cl,  F,  S),  species- sensitive  gas  chromatography 
(for  mercaptanes),  mass  spectrometry  (MS,  for  siloxanes 
and  chlorofluoro-hydrocarbons)  and  MS  combined  with 
GC  (for  complex  hydrocarbons).  Details  of  these  analyses 
are  outside  the  scope  of  this  paper  and  will  be  reported 
separately. 

The  model  was  defined  via  a  graphical  user  interface  in 
the  programme  “VAFI”  from  the  company  BEFSIM  (Fiege, 
Belgium).  VAFI  is  an  equation-based  data  reconciliation 
software,  including  an  extensive  thermodynamic  database, 


used  in  petrochemical,  chemical  and  power  plants.  Plant  data 
are  reconciled  and  computed  in  such  a  way  that  mass  and 
heat  balances  are  satisfied.  The  plant  under  study  is  imple¬ 
mented  in  a  process  flow  diagram,  using  streams  (gas,  liq¬ 
uid  or  both)  to  connect  equipment  units  (a  reactor,  a  heat 
exchanger,  a  pump,  etc.).  The  complete  model  constructed 
here  contained  a  total  of  192  variables  and  equations.  Typ¬ 
ical  computing  time  of  one  run,  i.e.  the  steady-state  output 
of  the  plant  for  a  given  set  of  input  parameters,  is  on  the 
order  of  seconds.  Output  is  stored  in  tags  and  given  as  (i) 
the  stream  compositions  at  any  position  in  the  process  flow 
scheme,  as  (ii)  the  heat  consumption  or  generation  of  each 
equipment  unit,  as  (iii)  temperature  inlets  and  outlets,  etc.  A 
graphical  MATFAB  interface,  programmed  in-house,  allows 
to  scan  the  variation  of  a  particular  parameter  (e.g.  the  CO2 
fraction  in  biogas,  the  biogas  flow  rate,  the  steam-to-carbon 
ratio,  the  excess  air  ratio  A  used  in  the  fuel  cell  stack,  etc.) 
and  record  the  output  of  every  computational  run.  Finally, 
output  is  also  represented  as  composite  curves  versus  heat 
flux  (kW). 

3.  Model 

The  process  flow  scheme  implemented  in  the  steady-state 
model  is  depicted  in  Fig.  1.  Streams  are  numbered  1-8.  Bio¬ 
gas  composition,  based  on  experimental  data,  is  precisely 
defined  in  stream  “1”,  which  is  preheated  and  mixed  with 
preheated  steam  (which  will  be  the  reference  case)  and/or 
air  for  conversion  to  synthesis  gas  in  the  reformer.  This  is 
carried  out  at  the  same  temperature  as  that  of  the  fuel  cell 
stack  (800  °C),  assuming  thermally  fully  integrated  reform¬ 
ing.  Cathode  air  is  preheated  to  an  inlet  temperature  of  min¬ 
imally  650  °C  (i.e.  a  maximum  of  150  K  below  the  stack 
temperature).  Inlet  air  is  50%  humid  and  contains  a  dry  frac¬ 
tion  of  79%  nitrogen.  Inlet  flows  are  circulated  by  blowers 
(1.1  bar,  i.e.  allowing  for  lOOmbar  pressure  drop  in  the  var¬ 
ious  system  units).  At  stack  temperature,  air  flow  is  divided 
into  a  pure  oxygen  stream  (“4a”)  deviated  to  the  anode — this 
separation  effectively  representing  the  solid  electrolyte — and 
the  excess  stream  “4b”,  which  combines  with  the  anode  out¬ 
let  (stream  “7”)  to  completely  convert  the  fuel  residue  to 
products  in  the  afterburner  zone,  the  adiabatic  temperature 
of  which  is  calculated.  Transferred  oxygen  (stream  “4a”)  and 
converted  fuel  (stream  “6”  minus  stream  “7”)  are  connected 
and  determined  by  a  fixed  value  for  the  fuel  utilisation  Uf , 
typically  set  to  80%.  All  input  streams  are  cold  fluxes  and 
require  heat  input  (the  g-arrows  in  Fig.  1).  The  exhaust  is 
the  only  hot  flux,  from  which  heat  is  recovered  to  preheat  all 
inputs.  The  stack  represents  a  hot  source  (heat  generation) 
and  the  reformer  a  cold  source  (heat  absorption),  both  at 
constant  temperature  (800  °C).  Electrical  power  is  removed 
from  the  stack,  consisting  of  SOFC  plates  20  cm  x  20  cm 
large  (361cm2  active  area  each)  connected  in  series.  The 
total  heat  balance  represents  the  thermal  efficiency  of  the 
system. 
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Fig.  1.  Process  flow  diagram  of  a  solid  oxide  fuel  cell  system  fed  with  sewage  biogas. 


3.1.  Electrochemistry 


The  electrochemical  model  for  fuel  conversion  in  the  stack 
was  elaborated  in  detail  in  a  previous  paper  [6].  Empirical 
expressions  for  cathode  and  anode  activation  polarisation 
loss  were  defined,  as  well  as  one  for  total  ohmic  losses 
(including  metal  interconnects  and  geometry  effects).  For 
the  present  study,  we  added  limitations  for  diffusion  at  anode 
and  cathode.  They  have  been  taken  into  account,  not  by 
using  diffusion  coefficients  applying  to  the  given  geometry 
[7],  but  by  the  basic  expressions  [8]: 


?7a  = 


Jlim  _ 


and 


0c  = 


Jlim  _ 


(la) 


(lb) 


for  anode  and  cathode  diffusion  overpotential  r]a  and  r]c,  re¬ 
spectively,  where  jnm  is  the  theoretical  current  density  for 
100%  fuel  or  oxygen  conversion  (X  denotes  the  cathode  air 
excess).  Whereas  this  description  oversimplifies  the  diffu¬ 
sion  process  in  fuel  cell  electrodes,  it  displays  the  merit  of 
correctly  “bending  down”  current-voltage  ( 1-V)  response  of 
the  fuel  cell  at  high  fuel  utilisation,  as  is  observed  experi¬ 
mentally.  In  particular  expression  (la)  for  the  anode  likely 
overestimates  the  diffusion  overpotential  loss — but  does  bear 
reference  to  the  fact  that  a  thick  support  electrode  of  limited 
porosity  is  assumed  to  be  used. 


All  voltage  losses  defined  in  the  model  are  simultane¬ 
ously  plotted  in  Fig.  2  as  a  function  of  current  density,  where 
in  order  to  plot  Eqs.  (la)  and  (lb),  which  require  knowl¬ 
edge  of  the  limiting  current  (i.e.  the  fuel  input  flow),  the 
latter  is  fixed  for  simplicity  at  0.5A/cm2.  An  appropriate 
current  density  regime,  where  long-term  stability  >40,000  h 
has  been  proven,  is  around  0. 2-0.4  A/cm2  [9,10].  High  cur¬ 
rent  density  on  smaller  cells  and  stacks  has  of  course  been 


Fig.  2.  Voltage  loss  equations  for  anode  and  cathode  used  in  the  electro¬ 
chemical  model  of  the  fuel  cell  stack,  based  upon  Ref.  [6]  and  Eqs.  (la) 
and  (lb).  R0hmic  =  0.13  Q  cm2.  To  illustrate  the  simple  diffusion  over¬ 
voltage  loss  (Eqs.  (la)  and  (lb)),  a  limiting  current  of  0.5  A/cm2  was 
taken  as  an  example  (X  =  3  for  the  airflow). 
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demonstrated,  also  for  considerable  operating  time  of  3000  h 
[11].  Nonetheless,  we  decided  to  maintain  the  conservative 
estimates  of  Fig.  2,  that  are  realistically  obtained  for  SOFC 
systems  today,  and  to  not  extrapolate  individual  best  data 
for  general  future  projections.  High  current  operation  also 
increases  heat  generation  on  the  stack,  implying  design  con¬ 
straints  for  cooling  (high  air  flow,  sufficient  reforming  en- 
dothermicity),  as  shown  further  on. 

3.2.  Carbon  deposition 


Likelihood  of  sooth  formation  from  the  fuel  mixture  is 
taken  into  account  by  considering  the  following  three  reac¬ 
tions: 


methane  pyrolysis  CH4  C  +  2H2  (2a) 

Boudouard  disproportionation  2CO  C  +  CO2  (2b) 

and  reverse  gasification  CO  +  H2  C  +  H2O  (2c) 


for  which  the  equilibrium  constants,  at  given  temper¬ 
ature  and  pressure  are  calculated  from  the  thermodynamic 
properties  of  the  reaction  compounds.  These  are  then 
compared  against  observed  values,  K ^  calculated  from  the 
actual  molar  fractions  (xz  )  and  total  pressure  (Ptot)  for  a  given 
stream: 


^kin  _ 

^  Rev  Gas 


(  PK2O  \  po 

\PU2PCoJ 


VH2Q  \  / 
XU2xCO  /  \  Ftot  / 


(3c) 


In  a  third  operation,  after  logarithmisation,  observed  val¬ 
ues  are  then  compared  against  the  theoretical  ones.  The  re¬ 
sult  of  the  operation: 

'  ^kin  ' 


then  indicates  whether  carbon  deposition  occurs  (result  <  0) 
or  not  (result  >  0),  for  each  of  the  three  equations.  This  cal¬ 
culation  was  performed  for  stream  “6”  in  the  model.  Care 
was  taken  to  always  observe  sufficiently  positive  values  for 
the  three  operations,  whenever  changing  inlet  fuel  compo¬ 
sition. 

Fig.  3  plots  the  equilibrium  values  versus  tempera¬ 
ture  for  the  three  reactions  (2a)-(2c).  It  clearly  illustrates  the 
critical  temperature  range  from  ^560  to  660  °C,  where  the 
tendency  for  carbon  formation  is  positive  for  all  three  re¬ 
actions  simultaneously  (or  in  other  words,  the  standard  free 
enthalpy  of  all  three  reactions  is  negative  in  this  temperature 
range). 


Fig.  3.  Equilibrium  constants,  Kp,  for  the  three  considered  carbon  forma¬ 
tion  reactions  (Eqs.  (2a)-(2c)),  as  a  function  of  temperature. 


4.  Results  and  discussion 

4.1.  Bio  gas  composition  and  production 

The  level  of  the  main  constituents  CH4,  CO2,  O2,  N2 
and  H2O  (vol.%)  as  well  as  of  H2S  (ppmv),  of  digestion 
gas  from  two  local  sewage  treatment  plants  was  measured 
regularly  during  several  months.  Oxygen  and  nitrogen  were 
not  separated  on  the  GC  column  and  are  therefore  taken 
as  a  sum  value  (“air  peak”).  Water  vapour  content  was  not 
measured  but  taken  as  the  difference  between  100%  and 
the  sum  of  the  signals  for  CH4,  CO2  and  air,  these  latter 
three  having  been  precisely  calibrated  on  known  premixed 
gases.  Table  1  summarises  the  data.  All  entries  represent 


Table  1 


Measured  biogas  compositions  of  two  sewage  plants  (~100kWei)  over 
the  period  July-November  2002 


%  ch4 

%  co2 

%  (N2  +  o2) 

%  h2o 

ppm  H2S 

Site  1 

63.33 

33.74 

1.14 

1.79 

1.52 

63.59 

33.13 

0.28 

3.00 

0.40 

62.16 

35.48 

0.87 

1.49 

1.21 

63.03 

33.85 

1.07 

2.05 

2.28 

63.55 

34.25 

0.75 

1.45 

3.77 

61.60 

36.37 

0.92 

1.11 

2.28 

63.28 

34.39 

1.30 

1.03 

2.42 

60.75 

36.79 

0.87 

1.59 

1.44 

64.42 

34.52 

1.06 

Dry 

— 

60.52 

35.78 

0.92 

2.78 

— 

Mean 

62.62 

34.83 

0.92 

1.81 

1.92 

Site  2 

60.10 

37.43 

0.21 

2.26 

0.66 

64.05 

35.08 

0.12 

0.75 

0.75 

62.47 

36.21 

0.10 

1.22 

2.15 

65.22 

32.10 

0.15 

2.53 

0.53 

63.67 

35.75 

0.57 

Dry 

— 

62.81 

35.78 

0.39 

1.02 

— 

Mean 

63.05 

35.39 

0.26 

1.56 

1.02 
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Table  2 

Production  data  of  an  existing  sewage  treatment  plant  averaged  over  the 
period  1997-2002 


Biogas 

Biogas  production  (m3/year)  378,439.5 

Used  in  the  130kW  engine  (%)  87.7 

Used  in  a  boiler  (%)  2.9 

Flared  off  %  9.3 

Equivalent  biogas  per  inhabitant  (1/day)  36.1 

Methane  equivalent  in  engine  (63%)  (m3/day)  571.8 

Electricity 

Yearly  engine  hours  (h/year)  5,613.3 

Average  daily  engine  hours  (h/day)  15.4 

Electricity  production  (MWhei/year)  463.0 

Average  load  (kWei)  82.7 

Methane  input  in  engine,  HHV  (J)  (MW  h/year)  2,187.0 

Methane  input  in  engine,  HHV  (W)  (kWei)  391.5 

Efficiency  HHV  (%)  21.2 

Site  consumption  (MWhei/year)  913.2 

Site  self-sufficiency  (%)  51.9 

Waste  water 

Waste  water  entry  (WW)  (m3/year)  4,591,667 

Fresh  sludge  (m3/year)  22,107 

Decanted  sludge  to  digester  (m3/year)  11,151 

Digested  sludge  out  (m3/year)  1,965 

Biogas  per  organic  dry  matter  (ODM)  (m3/t)  633.7 

Organic  charge  in  sludge  (%)  5.5 

FeCl3  consumption  (42%  solution)  (t/year)  284.7 

Fe  equivalent  (g/m3  WW)  8.6 


averages  of  minimum  20  consecutive  measurements  on  the 
same  sample  (total  analysis  time  ^3  h). 

From  Table  1,  we  derive  and  define  an  average  sewage 
biogas  composition  of  63%  CH4,  35%  CO2,  0.5%  air  and 
1.5%  H2O.  H2S,  which  is  noted  at  a  very  low  concentration 
of  1-2  ppm,  is  not  further  taken  into  consideration.  This 
biogas  would  not  require  sulphur  cleanup  (e.g.  by  active 
carbon  filters  or  ZnO  reactors)  before  admission  to  the  SOFC 
reformer  and  anode  catalysts. 

Actual  yearly  production  data  for  “Site  1”  (Table  1)  are 
given  in  Table  2.  Values  have  been  averaged  for  6  consecu¬ 
tive  years  (1997-2002). Yearly  biogas  production  is  an  aver¬ 
age  378,440  m3  (1036  m3/day),  of  which  ^88%  is  consumed 
by  the  gas  engine,  operating  for  about  ^15  h/day.  This  leads 
to  an  electricity  production  of  463  MW  h/year,  only  52%  of 
the  site  requirement  (913  MW  h/year).  Note  the  low  aver¬ 
age  efficiency  of  the  130kW  engine  (21%).  The  equivalent 
biogas  production  per  inhabitant,  connected  to  this  sewage 
network,  amounts  to  36 1/day.  At  63%  of  methane  content, 
this  figure  (8.3  m3  CHzt/year),  extrapolated  to  the  actual  Eu¬ 
ropean  Community  (EU-15,  375  mio  inhabitants),  leads  for 
this  region  to  an  ultimate  potential  of  2.67  Mtoe  (million 
tonnes  oil  equivalent),  corresponding  to  an  installed  elec¬ 
tricity  generation  potential  of  over  1  GW  (^10,000  sites 
of  ^lOOkWei).  Estimates  for  the  actually  exploited  sewage 
biogas  in  Europe  range  from  0.34  to  0.58  Mtoe  [12,13]. 

From  the  biogas  production  in  Table  2,  an  hourly  rate  of 
43  m3  is  derived,  which  is  used  in  the  SOFC  system  model  as 


input  value.  It  corresponds  to  a  higher  heating  value  (HHV) 
input  of  269  kW  (LHV  of  242  kW). 

No  data  on  the  heat  requirement  or  consumption  are  mea¬ 
sured  on  site.  Heat  is  recovered  from  the  engine  cooling  wa¬ 
ter  and  hot  exhaust  fumes,  to  provide  space  heating  and  to 
keep  the  digester  at  optimal  temperature  (35  °C).  This  suf¬ 
fices  to  cover  annual  thermal  needs,  with  a  large  surplus  in 
summer  (rejected)  and  only  a  very  slight  deficit  in  winter, 
during  which  a  total  of  around  1000  m3  of  natural  gas  (i.e. 
corresponding  to  <2  days  of  biogas  production)  is  burned 
in  the  boiler  for  supplemental  heating. 

Note  finally  the  important  consumption  in  iron  chloride 
solution,  added  to  the  incoming  waste  water  for  cold  desul- 
phurisation,  and  leading  to  the  low  sulphur  content  in  the 
biogas  product. 

4.2.  Determination  of  model  input  parameters 
4.2.1.  Reforming  conditions 

Steam  has  to  be  added  to  biogas  (in  this  case  the 
methane: carbon  dioxide  mixture  in  proportion  of  63:35) 
in  order  to  avoid  carbon  formation.  The  minimum  steam 
quantity  required,  at  a  given  temperature,  can  be  determined 
from  thermodynamic  equilibrium  concentration  calculation, 
and  is  plotted  in  Fig.  4  as  the  steam-to-carbon  ratio  (S/C, 
mole  H2O  versus  mole  C)  in  the  inlet  mixture.  At  800  °C, 
the  minimal  S/C  value  is  determined  to  0.37.  To  introduce 
a  safety  factor,  S/C  is  chosen  to  0.5  for  the  subsequent 
model  calculations,  effectively  corresponding  to  a  AT  of 
— 70K.  This  means  that,  in  practice,  the  endothermal  re¬ 
former  could  be  allowed  to  cool  from  800  to  730  °C  (during 
transients,  or  because  of  inertia),  without  running  into  ther¬ 
modynamic  risk  of  carbon  deposition.  Fig.  4  shows  that  for 
S/C  =  1.3,  carbon  formation  is  theoretically  excluded  at 
any  temperature. 


Fig.  4.  Minimal  temperature  and  corresponding  required  ratios  of  steam- 
to-carbon  (S/C)  and  oxygen-to -methane  (O2/CH4)  above  which  no  carbon 
deposition  takes  place  thermodynamically. 
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Also  displayed  in  the  graph  is  an  equivalent  condition 
when  using  partial  oxidation  as  fuel  processing  step  [6]. 
The  minimal  oxygen-to-methane  ratio  necessary  to  avoid 
sooth  building  from  the  63:35  CH4:C02  mixture  at  800  °C 
is  0.27.  As  partial  oxidation  is  fast  and  exothermic,  and  as 
we  cannot  allow  the  O2/CH4  ratio  to  reach  high  values  (to 
keep  the  electrical  efficiency  at  a  respectable  level),  O2/CH4 
is  limited  to  0.3.  This  choice,  corresponding  to  a  AT  of 
—20  K,  will  also  be  applied  as  syngas  formation  condition 
for  the  SOFC  system  model  (see  Section  4.5),  to  compare 
with  steam  reforming. 

4.2.2.  Reference  case 

Using  the  condition  S/C  =  0.5  at  the  fuel  inlet — in  the  fol¬ 
lowing  denominated  as  “reference  case” — thermodynamic 
equilibrium  concentrations  as  a  function  of  temperature 
(such  as  those  corresponding  to  the  reformer  outlet  stream, 
at  800  °C)  are  given  in  Fig.  5.  The  input  fuel  is  now  defined 
by  the  generic  formula  C0.98H3 .52O1.2N0.004S0.  At  each 
temperature  in  the  graph  this  atomic  balance  can  be  veri¬ 
fied.  At  the  reforming  temperature  of  800  °C,  conversion  to 
syngas  is  fairly  complete  and  the  H2/CO  ratio  around  1.8, 
in  between  that  of  pure  steam  reforming  (H2/CO  =  3)  and 
that  of  dry  reforming  (H2/CO  =  1)  since  we  effectively 
deal  with  mixed  reforming  of  CH4  with  both  H2O  and  CO2. 
This  is  simply  derived  from  the  nominal  reactions,  with 
S/C  =  0.5  and  a  CO2  fraction  of  35%: 

0.5 (CH4  +  H20  3H2  +  CO) 

0.35(CH4  +  C02  4^  2H2  +  2CO)  (  j 

adding  up  to: 

0.85CH4  +  0.5H2O  +  0.35CO2  ^  2.2H2  +  1.2CO  (4b) 


and  normalising  to  the  input  of  0.63  mole  methane  per  mole 
of  biogas: 

0.63CH4  +  0.37H2O  +  0.26CO2  1.63H2  +  0.89CO. 

(4c) 

The  last  equation  proves  the  H2/CO  ratio  of  1.63/0.89  = 
1.83  (Fig.  5)  and  also  that  0.12  mole  of  initial  steam  and 
0.09  mole  of  initial  carbon  dioxide  remain  unreacted,  as 
evident  in  the  figure.  Concentration  variations  in  FU,  H2O, 
CO  and  CO2  above  800  °C  (with  vanishing  CH4  and  absence 
of  carbon)  are  then  only  governed  by  the  shift  reaction. 

The  dotted  line  indicates  the  open  circuit  voltage  (OCV) 
that  would  be  observed  for  a  perfectly  sealed  cell  with  the 
anode  inlet  gas  composition  defined  in  Fig.  5  at  each  corre¬ 
sponding  temperature  with  respect  to  air  ( po2  =0.21  atm) 
at  the  cathode.  OCV  at  800  °C  is  then  1.05  V. 

The  stoichiometric  oxygen  amount  needed  to  burn  all  fuel 
is  given  by: 

no2  nu  no 

— -  =nc-\-  — - — 

nf  4  2 

=  0.98  +  0.88  -  0.6  =  1.26  ( uf  =  100%),  (5) 

obviously  leading  to  2  mole  equivalent,  of  oxygen  per  mole 
of  methane  (2  x  0.63  =  1.26)  as  essentially  no  free  oxygen 
is  present  in  the  fuel  (0.1%,  which  is  neglected). 

Adding  oxygen  to  the  reformed  fuel  (stream  “4a”  to 
stream  “6”),  anodic  oxidation  proceeds  and  the  outlet  com¬ 
position  is  calculated  as  a  function  of  fuel  utilisation,  uf , 
plotted  in  Fig.  6,  showing  the  consumption  of  H2  and  CO 
fuels  and  formation  of  H2O  and  CO2  products.  Also  shown 
is  the  corresponding  OCV,  dropping  to  0.87  V  at  the  outlet 
for  80%  fuel  utilisation.  As  the  model  computes  the  stack 
electrical  output  using  the  electrochemical  losses  defined 
in  Fig.  2,  that  are  subtracted  from  the  free  enthalpy  of  the 


Fig.  5.  Thermodynamic  equilibrium  concentrations,  as  a  function  of  tem¬ 
perature,  for  the  fuel  mixture  consisting  of  1  mole  biogas  (63%  CH4, 
35%  CO2,  1.5%  H2O,  0.5%  N2)  with  0.475  mole  water,  corresponding 
to  an  S/C  ratio  of  0.5.  Molar  concentrations  normalised  to  0.98  mole  of 
carbon  in  the  mixture. 


fuel  utilisation  u  F  (-) 


Fig.  6.  Thermodynamic  equilibrium  concentrations,  as  a  function  of  fuel 
conversion  (%),  i.e.  as  a  function  of  0-90%  oxygen  addition,  for  the 
product  mixture  of  1  mole  biogas  (63%  CH4,  35%  CO2,  1.5%  H2O,  0.5  % 
N2)  and  0.475  mole  water  at  800  °C  (see  Fig.  5).  Molar  concentrations 
normalised  to  0.98  mole  of  carbon. 
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reaction,  in  itself  determined  by  anode  plus  cathode  outlets 
minus  inputs  (OCV  at  the  inlet  1.05  V,  at  the  outlet  0.87  V), 
an  average  OCV  value  is  to  be  considered,  i.e.  0.96  V. 


4.2.3.  Air  excess  ratio 

A  multiple  quantity  of  the  stoichiometric  oxygen  amount, 
indicated  by  the  factor  A,  is  flown  at  the  cathode  for  stack 
cooling.  The  overall  fuel  conversion  reaction  is  therefore 
written: 


C0.98H3.52O1.2N0.005  +  1.26A(02  +  3.773N2) 

=>  0.98CO2  +  1.76H20  +  4.754N2  +  1.26(A  -  1)02. 

(6) 


When  A  =  1,  the  water  vapour  fraction  in  the  outlet  is 
maximal  and  equal  to  23.5%.  With  an  empirical  formula  for 
the  saturation  pressure: 


Psat(Pa)  =  exp 


23.57 


4041.6 

rsat(°C)  +  235.6 


we  obtain  64  °C  for  the  temperature  at  which  water  vapour 
will  start  to  condense  from  the  mixture.  Usually,  A  is  higher, 
so  Tsat  will  be  lower  still.  Nevertheless,  outlet  temperature 
for  the  cold  exhaust  fumes  (Fig.  1,  stream  “8b”)  was  fixed 
to  a  realistic  80  °C,  without  recovering  condensation  heat. 

Given  the  biogas  inlet  flow  (43  m3/h),  the  fuel  utilisation 
(0.8)  and  the  electrochemical  losses  (Fig.  2),  the  necessary 
number  of  cells  for  the  stack  can  be  estimated.  Before,  two 
further  assumptions  are  made.  First,  A  is  for  the  moment  pos¬ 
tulated  at  the  factor  3,  a  common  value  for  steam-reformed 
methane  driven  stacks  [14].  Second,  a  minimal  cell  operating 
voltage  of  0.65  V  is  postulated.  We  indeed  observe  experi¬ 
mentally  that  cell  degradation  accelerates  below  such  operat¬ 
ing  voltage,  which  can  be  related  to  the  reduction-oxidation 
potential  of  the  nickel  anode  catalyst  [15]. 

Following  the  previous  arguments,  a  cell  loss  of  around 
3 10  mV  can  therefore  be  tolerated  (0.96  V  OCV,  —0.65  V 
operating  voltage).  Combined  diffusion  losses,  from 
Eqs.  (2a)  to  (2c)  with  Uf  =  0.8  and  A  =  3,  can  be  de¬ 
termined  to  89  mV.  Hence,  activation  and  ohmic  loss  can 
total  approximately  220  mV.  From  Fig.  2,  we  then  obtain  a 
corresponding  current  density  of  0.24  A/cm2,  or  86.4  A  on 
the  stack  (361  cm2  active  area  per  cell). 

43m3/h  of  biogas  (63%  methane)  correspond  to  an  in¬ 
put  flow  of  0.5  mol/s.  Stoichiometric  oxygen  implies  1.26 
times  this  flow,  but  Uf  is  fixed  to  0.8.  Therefore,  the  oxygen 
flow  through  the  electrolyte  membrane  (stream  “4a”)  is  0.51 
mole/s.  Converting  to  current  (*4*F),  this  equals  196,560  A 
for  one  hypothetical  cell.  With  the  stack  current  determined 
to  86.4  A,  the  ratio  of  both  values  gives  the  required  cell 
number,  i.e.  ^2400. 

This  value  is  very  large  and  due  to  the  purposefully  cho¬ 
sen  conservative  loss  characteristics  (Fig.  2),  equivalent  to  an 
area  specific  resistance  (asr)  of  around  1.2  £2  cm2.  Our  own 
ASE  cells  presently  show  an  asr  of  approximately  0.6  Q  cm2 
[16],  so  that  half  the  amount  of  cells  (i.e.  1200)  could  even¬ 


AT  on  stack  (T  -  T.  )  (K) 

v  out  in'  v  ' 


Fig.  7.  Heat  sinks  (kW)  provided  by  endothermal  reforming  and  circulating 
gases,  as  a  function  of  the  air  flow  (A.)  and  of  the  thermal  gradient  on 
the  stack  (AT),  to  evacuate  the  heat  generated  in  the  fuel  cell  stack. 


tually  suffice  to  construct  the  100kWei  sized  stack,  with 
a  obvious  enormous  impact  on  cost.  For  comparison,  the 
demonstrated  100kWei  tubular  SOFC  stack  [4]  used  1152 
tubes  of  ^1000  cm2  active  area  each. 

Excess  air  A  had  been  postulated  to  the  factor  3  above.  In 
the  following,  justification  for  this  choice  is  given  through 
consideration  of  the  stack  thermal  balance.  Heat  developed 
on  the  stack  (gstack)  because  of  the  exothermal  reactions 
has  to  be  evacuated  to  avoid  thermal  runaway.  This  can 
be  accomplished  by  the  endothermal  reformer,  in  thermal 
conctact  with  the  stack,  on  the  one  hand,  and  on  the  other 
hand  by  the  gases  transiting  the  stack,  especially  the  cathode 
air  excess  flow.  Fig.  7  gives  the  evacuated  heat  (in  kW)  for 
the  reference  case  considered  (43m3/h  biogas  input  flow 
with  S/C  =  0.5).  In  particular,  heat  flux  required  to  heat 
anode  and  cathode  (depending  on  A)  gases  over  a  certain 
temperature  difference  AT  are  plotted.  It  appears,  once  the 
fuel  input  flow  and  the  reforming  mixture  are  fixed,  that 
the  air  excess  and  the  allowable  AT  are  the  only  control 
variables  that  can  regulate  heat  removal  from  the  stack,  the 
small  anode  gas  flow  playing  a  minor  role.  A  maximum 
AT  of  200 K  is  admissible  [17,18].  This  may  be  brought 
in  relation  to  thermal  stress  of  the  ceramic  materials  that 
compose  the  stack.  Stress  accumulation  in  a  material  due  to 
a  thermal  gradient  may  in  a  simple  form  be  expressed  as 
[19]: 

a(MPa)  =  aATEy,  (8) 

with  a  the  thermal  expansion  coefficient  in  m/m  K  and  E\ 
the  materials  Young  modulus  in  Pa.  Using  order  of  magni¬ 
tude  values  of  a  =  10~5  K_1  and  Ey  =  1011  Pa  for  the  ce¬ 
ramics,  a  AT  of  200  K  leads  to  thermal  stress  of  200  MPa, 
close  to  the  modulus  of  rupture  of,  for  example,  the  elec¬ 
trolyte  material  at  high  temperature  [20].  Similarly,  an  tol¬ 
erable  diagonal  gradient  on  the  proposed  planar  ASE  cells 
has  been  reported  to  0.7K/mm  [21],  hence  approximately 
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Fig.  8.  Generated  electrical  power  and  heat  (kW)  by  the  fuel  cell  stack,  as  a 
function  of  air  flow  (A).  Also  shown  as  a  function  of  A  is  the  cooling  power 
(sum  curve  from  Fig.  7  with  AT  =  200 K),  as  well  as  the  hypothetical 
adiabatic  temperature  in  the  postcombustion  zone.  The  reference  case  of 
steam-reformed  biogas  (S/C  =  0.5)  is  taken,  rstack  =  ^reformer  =  800  °C. 


200  K  on  a  20  cm  x  20  cm  square  cell.  For  this  reason,  air 
at  the  cathode  inlet  in  the  system  model  may  be  as  low  as 
650  °C  as  minimal  temperature,  and  the  outlet  (streams  “4b” 
and  “7”)  up  to  850  °C  at  most,  in  order  to  limit  materials 
deterioration.  Especially,  the  cheap  metal  interconnects  used 
for  cell  stacking  (FeCr  alloys)  cannot  tolerate  temperatures 
above  850  °C  in  an  atmosphere  with  high  steam  levels  for  a 
long  time  [22]. 

Fig.  8  displays  the  balance  between  heat  generation  on 
and  cooling  of  the  SOFC  stack.  For  the  chosen  input  condi¬ 
tions  (discussed  higher  and  summarised  further  in  Table  3), 
electrical  and  thermal  power  developed  on  the  stack  are 
131  kWei  and  143kWth,  depending  somewhat  on  the  circu¬ 
lating  air  excess.  The  curve  denoted  “stack  cooling  power” 
adds  up  the  evacuated  heat  fluxes  (Fig.  7)  by  the  reformer 
(78  kW)  and  the  cathode  air  for  a  AT  equal  to  200  K  (<2air  = 
A.  x  21  kW).  Stack  heating  and  cooling  are  exactly  balanced 
at  X  =  3. 

The  preceding  argument  of  course  is  a  rather  artificial 
construction.  In  physical  reality,  the  stack  may  experience 
further  heating  from  an  afterburner  zone  in  proximity  and 
cooling  because  of  insulation  losses.  However,  the  simple 
representation  used  here  illustrates  the  general  argument  of 
fuel  cell  stack  thermal  balancing  by  easily  traceable  numbers 
in  a  simplified  system  component  configuration. 

4.3.  Model  output:  steam-reformed  bio  gas 

Having  in  the  previous  sections  defined,  and  motivated 
the  selection  of,  all  input  parameters,  the  steady- state  model 
output  for  the  reference  case  (steam-reformed  biogas,  S/C  = 
0.5)  is  discussed  in  the  following.  Table  3  lists  all  relevant  re¬ 
sulting  numbers,  and  Table  4  gives  the  stream  compositions 
at  different  positions  in  the  process  flow  scheme,  whereas 


Table  3 


Input  and  output  data  of  the  process  flow  model  for  the  reference  case 
of  steam-reformed  biogas  (S/C  =  0.5) 


Biogas  inlet  flow  (Nm3/h) 

43 

Methane  inlet  flow  (Nm3/h) 

27.1 

Fuel  input  power  (kW) 

269.3 

Stack  inlet  temperature,  air  (°C) 

650 

Reformer  temperature  (°C) 

800 

Stack  temperature  (°C) 

800 

Cell  area  (cm2) 

361 

Cell  number 

2,400 

Air  excess  A 

3 

Steam-to-carbon  ratio  S/C 

0.5 

Stack  voltage  (V) 

1,553.9 

Cell  voltage 

0.648 

Current  (A) 

86.37 

Current  density  (A/cm2) 

0.2393 

7  a,  act  00 

0.0965 

7a,diff  (V) 

0.0744 

7c, act  (V) 

0.0884 

7c,  diff  00 

0.0143 

Gross  electric  power  (kWei) 

134.21 

Air  compressor  power  (kWei) 

-2.986 

Fuel  compressor  power  (kWei) 

-0.170 

Net  electric  power  (kWei) 

131.06 

Gross  electrical  efficiency  (%) 

49.83 

Net  electrical  efficiency  (%) 

48.66 

Thermal  efficiency  (%) 

39.57 

Stack  generated  heat  (kW) 

142.95 

Reformer  cooling  power  (kW) 

-78.24 

Stack  cooling  air  (AT  =  150  K)  (kW) 

-48.05 

Hot  exhaust  (kW) 

317.96 

Main  air  preheater  (kW) 

-187.24 

Fuel  preheater  (kW) 

-21.72 

Economiser  (kW) 

-1.70 

Evaporator  (kW) 

-10.49 

Steam  preheater  (kW) 

-6.88 

Cold  exhaust  loss  (kW) 

20.45 

Fig.  9  plots  the  stack  I-V  behaviour  at  design  fuel  input 
(43  m3/h). 

Table  3  shows  that  1 34  kWei  are  obtained  from  the  stack, 
leading  after  subtraction  of  3  kW  required  to  circulate  the 
gases  (850m3/h  at  1.1  bar  inlet)  a  net  dc  generation  of 
131  kWei,  at  an  HHV  efficiency  of  48.7%.  The  gas  engine 
installed  on  site  is  precisely  of  130kWei  size.  Generated 
heat  in  the  stack  is  143  kWth,  balanced  by  the  reformer  and 
the  cathode  air  as  discussed  (Fig.  8).  Available  heat  from  the 
exhaust  when  cooling  from  the  afterburner  outlet,  986  °C, 
to  the  released  fumes  at  80  °C,  is  318kWth,  sufficient  to 
preheat  air  (187kWth)  to  stack  inlet  temperature  (650  °C), 
fuel  (22kWth)  and  steam  from  reforming  water  (19kWth). 
The  total  heat  balance  amounts  to  +  106kWth  or  almost 
40%  thermal  efficiency  (269  kW  methane  input). 

Table  4  shows  the  reformer  outlet  is  90%  syngas  at  nearly 
double  the  fuel  inlet  flow  (116m3/h  versus  64m3/h),  of 
which  80%  is  converted,  or  1 16  x  0.9  x  0.8  =  86.6  m3/h,  re¬ 
quiring  just  half  of  this  in  oxygen  flow  through  the  electrolyte 
membrane  (stream  “4a”,  43.3  m3/h).  This  stream,  12102/s 
through  2400  cells  x  0.0361  m2  =  87  m2  of  active  area,  of 
course  translates  to  0.24  A/cm2  current  density. 
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Table  4 

Stream  compositions  at  various  positions  in  the  process  flow,  for  the  reference  case  of  steam-reformed  biogas  (S/C  =  0.5),  given  in  total  flow  (N  m3/h) 
and  in  mol%  per  stream 


Position  (mol%) 

Cathode 
flow  (m3/h) 

%  02 

%  n2 

%  h2o 

Fuel  flow 
(m3/h) 

%  ch4 

%  C02 

%  h2o 

%  n2 

%  h2 

%  CO 

Inlet 

780.5 

21 

78.1 

0.9 

63.9 

42.7 

23.7 

33.2 

0.4 

— 

— 

Reformer  outlet 

116.2 

0.78 

3.24 

5.7 

0.2 

57.9 

32.2 

O2  transfer 

43.3 

100 

Stack  outlet 

737.2 

16.2 

83 

0.9 

118 

0 

29.2 

52.2 

0.2 

11.9 

6.5 

Postcombustion  outlet 

844 

12.8 

72.4 

5 

9.8 

Fig.  9  compares  the  conservative  I-V  output  from  the 
model  with  that  of  recent  results  obtained  in-house  on  an¬ 
ode  supported  electrolyte  cells,  from  an  8 -cell  stack  test  of 
50  cm2  active  area  each,  delivering  140  Wei  below  800  °C 
(0.32  W/cm2)  from  H2  fuel,  however,  for  lower  fuel  conver¬ 
sion  (50%).  The  comparison  illustrates  the  realistic  perfor¬ 
mance  estimate  for  the  model  stack.  Efficiency  reaches  49% 
and  is  maximal  at  85%  fuel  utilisation,  after  which  anode 
diffusion  loss  causes  the  operating  voltage  to  drop  steeply. 

From  another  point  of  view,  input  methane  of  269  kW 
is  upgraded  to  347  kW  in  the  reforming  process  (78  kW 
input),  of  which  80%  (i.e.  277  kW)  is  converted  on  the 
stack,  namely,  38%  to  electricity  (134kW)  and  42%  to  heat 
(143  kW). 

Fig.  9  is  replotted  in  Fig.  10,  converting  current  density  to 
fuel  utilisation,  but  focussing  now  on  the  thermal  balance  of 
the  system,  indicated  in  kW.  Electrical  power  and  exploitable 
excess  heat  (“thermal  power”)  always  add  up  to  238  kW 
(88.5%  cogeneration  efficiency),  with  the  electricity  fraction 
dominating  the  heat  fraction  once  above  65%  fuel  utilisation. 
While  electricity  production  flattens  out  (126-132  kWei)  be¬ 
tween  70  and  90%  fuel  conversion,  heat  generation  in  this 
regime  sharply  increases  (115-1 80  kW),  requiring  an  adap¬ 


current  density  /  Acm 


Fig.  9.  Current  density-cell  voltage-power  density  output  for  the  model 
cell  stack  (800  °C),  considering  the  reference  case  of  biogas  feed  with 
S/C  =  0.5.  Electrical  efficiency  is  seen  to  raise  to  49%,  at  80%  fuel 
utilisation  and  0.15  W/cm2  power  density.  Also  shown  is  the  experimental 
cell  output  of  an  in-house  developed  stack  (140  Wei,  800  °C,  Wf,  50%), 
illustrating  adequate  realistic  potential. 


tation  in  cathode  air  flow  if  materials  exposure  to  excessive 
local  temperature  is  to  be  avoided.  Heat  evacuation  by  re¬ 
forming  and  cathode  air  (A.  =  3,  AT  =  200  K),  totalling 
142  kW,  is  also  displayed. 

4.4.  Composite  curves 

The  following  three  graphs  (Figs.  11-13)  plot  the  com¬ 
posite  curves  for  the  reference  case  (S/C  =  0.5)  of  the  bio¬ 
gas  SOFC  system  model. 

In  Fig.  11,  the  upper  curve  or  hot  composite  is  the  sum 
of  all  hot  fluxes  and  sources  in  the  system,  sorted  by  tem¬ 
perature  level:  Gexhaust  from  the  postcombustion  tempera¬ 
ture  (986  °C  or  1259  K)  down  to  800  °C,  Gstack  as  source 
at  800  °C,  and  Gexhaust  again  from  800  °C  down  to  the  final 
temperature  (15  °C)  of  rejection  to  the  environment  (which 
includes  “Gloss’’  from  80  to  15  °C,  see  Fig.  1).  The  bottom 
part  of  the  hot  composite  fixes  the  scale  of  exchanged  heat 
(v-axis,  Q  in  kW)  at  zero.  In  total,  481.5  kW  is  available,  as 
verified  from  Table  3  (Gstack  =  143  kW,  Gexhaust  =  318  kW, 
Gloss  =  20.5  kW).  From  the  upper  part  of  the  hot  compos¬ 
ite  (1259  K,  481.5  kW),  the  heat  requirements  of  the  system 
are  satisfied  “going  backwards”,  i.e.  the  cold  composite  or 


fuel  utilisation  u  F  (-) 


Fig.  10.  Heat  balance  of  the  fuel  cell  stack,  as  a  function  of  fuel  conver¬ 
sion.  Generated  electrical  power  and  useful  thermal  output  power  balance 
out  to  a  constant  total  output  of  approximately  238  kW  (88.5%  HHV  co¬ 
generation  efficiency),  but  generated  stack  heat  strongly  increases  at  high 
fuel  conversion,  requiring  appropriate  cooling  (by  endothermal  reforming 
and  cathode  excess  air,  here  k  =  3). 


136 


J.  Van  herle  et  al.  / Journal  of  Power  Sources  131  (2004)  127-141 


Composite  curves  of  the  system 


Fig.  11.  Hot  and  cold  composite  curves  of  the  SOFC  biogas  system  (S/C  =  0.5).  Plateaus  of  stack  heat  generation,  endothermal  reforming  and  steam 
generation  are  easily  recognised.  Available  excess  heat  (final  difference  between  the  two  curves)  is  126  kW. 


the  sum  of  all  cold  fluxes  and  utilities  in  the  system  (lower 
curve  in  Fig.  1 1)  is  traced  back  towards  the  left  on  the  v-axis. 
These  are:  ^reformer  (78  kW)  as  cooling  utility  (at  800  °C, 
but  shown  displaced  in  the  figure,  at  a  level  of  790  °C,  for 
clarity),  preheating  of  air,  steam  and  fuel  between  100  and 
800  °C,  evaporation  of  water  to  reforming  steam  (10.5  kW), 
and  preheating  of  all  fluxes  between  inlet  temperature  (15  or 
25  °C)  and  100  °C.  The  needs,  as  verified  from  Table  3  (i.e. 
all  negative  heat  fluxes  Q ),  total  355  kW.  The  difference  be¬ 
tween  both  hot  and  cold  composite  streams  equals  126.5  kW, 
visible  on  the  v-axis  between  origin  and  low  temperature 
end-point  of  the  cold  stream.  Subtracting  the  final  loss  to  the 
environment  (20.5  kW),  the  available  cogeneration  heat  flux 
(106  kW,  thermal  efficiency  39.7%)  is  obtained.  This  heat  is 
available  at  a  temperature  level  at  least  of  the  hot  composite 
flux,  at  which  all  cold  flux  preheating  is  satisfied,  i.e.  650  K 
in  the  figure. 

Fig.  12  shows  the  grand  composite  curve  of  the  system 
versus  heating  power  (kW),  i.e.  the  difference  between  hot 
and  cold  composite  fluxes  (Fig.  11)  or  net  heat  flux  between 
each  two  corresponding  temperature  levels.  The  origin  of 
the  v-axis  is  now  placed  at  the  highest  temperature  level 
(1259  K),  starting  with  the  hot  exhaust  from  the  afterburner, 
and  plotting  net  heat  generation  towards  the  right  and  net 
heat  absorption  towards  the  left.  Again  the  three  isothermal 


plateaus  corresponding  to  gstack  (generation  of  143  kW,  at 
800  °C),  ^reformer  (absorption  of  78  kW,  drawn  at  790  °C, 
for  clarity)  and  ^evaporator  (absorption  of  lOkW,  at  100  °C) 
are  found.  The  grand  composite  shows  that  the  available  net 
heat  flux  of  126  kW  (cogeneration  heat),  is  actually  available 
at  the  temperature  of  800  °C,  i.e.  at  very  high  exergy  value. 
Therefore,  in  reality,  the  overall  electrical  efficiency  of  the 
system  could  still  be  enhanced  by  pressurising  the  flows  and 
adding  an  expansion  microturbine  to  the  exhaust. 

Fig.  13  finally  shows  another  alternative  to  Fig.  11,  now 
plotted  with  the  exergy  equivalent  on  the  y-axis  instead  of  ab¬ 
solute  temperature.  The  factor  (1-288  K IT  (K))  corresponds 
to  the  Carnot  efficiency  factor  for  converting  heat,  between 
the  two  temperature  levels  T  and  To  (288  K,  the  environ¬ 
ment),  into  work.  The  zone  in  between  both  exergy  compos¬ 
ite  curves  indicates  all  exergy  loss  occurring  in  internal  heat 
exchange,  if  the  present  configuration  of  the  system  were  to 
be  maintained  as  such.  In  other  words,  the  diagram  shows 
the  potential  for  improved  heat  exchange  within  the  sys¬ 
tem,  that  could  be  achieved  when  approaching  both  curves 
as  closely  as  possible.  This  would  require  further  close  in¬ 
vestigation  of  the  physical  layout  of  the  various  heat  ex¬ 
changes  taking  place,  attributing  appropriate  exchange  areas 
and  transfer  coefficients  for  each  case.  This  will  be  reported 
in  a  forthcoming  study. 
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Fig.  12.  Grand  composite  curve,  or  difference  between  hot  and  cold  composite  curve  of  Fig.  11  of  the  SOFC  biogas  system,  showing  available  excess 
heat  at  high  temperature  level  (800  °C). 


Composite  curves  of  the  system 


/V  _ 

4.5.  Operating  parameter  sensitivity 
4.5.1.  Variable  load 

SOFC  stack  behaviour  with  variable  biogas  input  flow  is 
discussed  in  the  following.  At  the  outset,  a  constant  flow 


of  43  m3/h  (design  point)  was  postulated,  in  reality  this  in¬ 
put  flow  may  fluctuate  according  to  demand  for  electricity 
or  heat  or  else  according  to  biogas  production  availability. 
Results  for  a  variable  load,  i.e.  for  changing  biogas  input 
flow  between  20  and  200%  of  the  design  value,  are  given 


Fig.  13.  Exergy  hot  (upper)  and  cold  (lower)  composite  curves  of  the  SOFC  biogas  system. 
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Fig.  14.  SOFC  stack  behaviour  (efficiency,  I-V)  as  a  function  of  changing 
biogas  input  flow  (20-200%  of  the  design  point  43m3/h).  Electrical 
efficiency  raises  at  partial  load  and  drops  at  higher  load. 

in  Figs.  14  and  15.  Fig.  14  shows  quantitatively  how  elec¬ 
trical  efficiency  drops  (between  61  and  37%)  with  increas¬ 
ing  load,  due  to  higher  current  density  on  the  cells  (up  to 
0.48  A/cm2)  at  lower  operating  voltage  (down  to  0.51  V  per 
cell).  In  effect,  we  move  up  and  down  the  stack  I-V  curve 
when  changing  the  fuel  flow,  while  keeping  the  conversion 
(80%)  constant.  Fig.  15  shows  that  the  electrical  output  be¬ 
tween  20  and  200%  fuel  load  varies  from  30  to  205  kW 
(design  point  =  131  kW).  At  the  same  time,  even  though 
air  excess  A  has  been  kept  constant  (A  =  3),  i.e  cathode 
air  increases  concomitantly  when  fuel  input  increases,  the 
thermal  gradient  on  the  stack  also  increases,  from  the  maxi¬ 
mum  tolerable  200  K  between  outlet  and  inlet  at  design  point 
(131  kW,  100%  load)  to,  for  example,  250  K  when  running 
at  150%  load  (65  m3  biogas/h).  In  other  words,  A  should  be 
further  increased  at  higher  load  and  may  be  decreased  at 
partial  load. 


Fig.  15.  SOFC  stack  output  and  thermal  gradient  as  a  function  of  changing 
biogas  input  flow  (20-200%  of  the  design  point  43  m3/h).  Power  increase 
is  not  linear  and  thermal  stress  increases  at  higher  load. 


4.5.2.  Fuel  composition 

The  final  graphs  illustrate  the  effect  of  fuel  composi¬ 
tion  variation  at  constant  fuel  inlet  flow  (43  m3/h),  either  by 
changing  the  reforming  steam  level  (S/C  ratio)  or  by  fluc¬ 
tuations  in  the  biogas  CO2  content.  While  the  composition 
of  sewage  biogas  is  stable  (Table  1),  making  it  one  of  its 
attractive  features,  the  composition  of  other  biogas  produc¬ 
tion  sites  may  fluctuate  strongly.  Fig.  16  shows  examples 
of  this,  the  first  (a)  was  measured  on  a  200  kW  solid  green 
waste  methanisation  plant  [23],  the  second  (b)  was  measured 
on  a  lOkW  farm  waste  digestion  installation  [23].  In  the 
first  case,  connected  engines  will  stop  functioning  once  the 
methane  level  drops  below  50%,  which  occurs  periodically. 
This  would  not  apply  to  a  fuel  cell.  In  the  second  case,  im¬ 
portant  transients  are  periodically  observed  (air  entry  due  to 
batch-refilling  of  the  digesters)  causing  a  simultaneous  drop 
in  methane  and  CO2  levels. 

Fig.  17  illustrates  that  an  increasing  CO2  level  in  the  bio¬ 
gas  stream  at  constant  flow  rate  obviously  decreases  power 
generation  on  the  SOFC  stack  (by  — 2kWei  per  1%  point  in¬ 
crease  in  CO2),  but  at  the  same  time  that  efficiency  slightly 
increases.  As  fuel  conversion  is  constant  (80%),  less  methane 


Fig.  16.  (a)  CH4  content  measured  in  biogas  from  a  200  kW  methanisation 
plant  of  solid  green  waste,  during  two  different  10-day  periods  (November 
2002,  February  2003).  (b)  CH4  and  CO2  content  measured  in  biogas  from 
a  lOkW  farm  waste  digester  during  a  1-year  test. 
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C02  content  (%)  in  biogas 


Fig.  17.  Variation  of  generated  stack  power  and  efficiency  with  the  CO2 
level  in  the  biogas  feed,  at  constant  flow  rate.  Power  varies  by  about 
2kWd  per  %  variation  in  CO2,  and  efficiency  slightly  drops  for  richer 
fuel  as  current  density  (constant  conversion,  80%)  then  increases. 

at  the  inlet  corresponds  to  smaller  current  density  on  the 
same  stack  and  hence,  moving  up  the  I-V  curve,  to  higher 
operating  voltage  and  efficiency. 

Fig.  18  shows  the  effect  of  further  steam  addition  to  the  re¬ 
former.  Doubling  the  steam  amount  (S/C  =  1)  of  the  design 
value  (S/C  =  0.5),  and  therefore  almost  excluding  carbon 
formation  anywhere  in  the  system  (Fig.  4)  lowers  the  elec¬ 
trical  output  from  131  to  127  kW  and  efficiency  (48.7-47%) 
by  only  a  small  amount.  While  current  density  remains  con¬ 
stant,  cell  voltage  does  slightly  drop  owing  to  higher  pu2o 
at  the  anode.  However,  thermal  efficiency  also  drops  (from 
39.7  to  36%),  due  to  lower  postcombustion  temperature  and 
higher  steam  preheating  requirement,  to  a  total  cogeneration 
HHV  efficiency  of  83%  (compared  to  89%  at  design  point, 
S/C  =  0.5). 


Fig.  18.  Variation  of  generated  stack  power  and  efficiency  with  S/C  ratio, 
both  showing  only  a  limited  drop  with  the  increase  in  added  steam  for 
reforming,  due  to  a  minor  lowering  in  cell  voltage  (higher  pu^o  at  the 
anode). 

For  final  evaluation,  steam-reformed  biogas  was  com¬ 
pared  to  a  number  of  other  input  fuels  or  reforming  con¬ 
ditions,  maintaining  the  same  system  layout  (Fig.  1)  and 
electrochemical  model.  These  are:  steam-reformed  pure 
methane,  partially  oxidised  biogas,  pure  hydrogen  and  pure 
methane.  The  latter  case  is  only  hypothetical  as  severe  car¬ 
bon  deposition  would  occur.  It  is  included  for  theoretical 
reference  purpose.  In  all  cases,  fuel  flow  was  adjusted  to 
give  identical  inputs  of  269.3  kW  everywhere.  Table  5  sum¬ 
marises  the  computed  outputs.  Steam-reformed  biogas  is 
essentially  identical  to  steam-reformed  methane;  the  CO2 
reforming  compound  inherent  in  biogas  increases  endother- 
mal  reformer  cooling  (78  kW  versus  71  kW),  therefore, 
upgrades  the  input  fuel  to  a  higher  heating  value  at  the 
reformer  outlet,  but  this  advantage  is  rather  exactly  compen- 


Table  5 


Comparison  of  stack  output  for  different  fuel  processing  methods  and  different  fuel  feeds 


Parameter 

Steam-reformed  biogas 

Biogas  POX 

Steam-reformed  CH4 

h2 

ch4 

Steam-to-carbon  ratio 

0.5 

— 

1.15 

— 

— 

02-to-CH4  ratio 

— 

0.3 

— 

— 

— 

Cell  voltage  (V) 

0.648 

0.682 

0.654 

0.697 

0.765 

Current  density  (A/cm2) 

0.239 

0.203 

0.239 

0.198 

0.239 

Gross  electric  power  (kWei) 

134.2 

120.3 

135.4 

119.6 

154.6 

Air  compressor  power  (kWei) 

2.99 

4.32 

2.98 

4.94 

1.99 

Net  electric  power  (kWei) 

131.1 

115.65 

132.29 

114.3 

152.6 

Gross  electrical  efficiency  (%) 

49.84 

44.67 

50.35 

44.55 

57.52 

Net  electrical  efficiency  (%) 

48.66 

42.94 

49.20 

42.58 

56.74 

Thermal  efficiency  (%) 

39.58 

46.36 

36.98 

45.63 

38.40 

Stack  generated  heat  (kW) 

142.9 

117.9 

135.4 

101.2 

44.3 

Reformer  cooling  power  (kW) 

78.2 

29.4 

71.1 

— 

— 

Air  excess  k 

3 

4.5 

3 

6 

2 

Stack  cooling  air  (kW) 

64.1 

92.9 

64.0 

106.5 

42.7 

AT  (K) 

202 

191 

201 

190 

207 

Hot  exhaust  (kW) 

318.0 

411.0 

313.2 

437.1 

230.2 

Cold  exhaust  loss  (kW) 

20.4 

29.1 

20.2 

31.8 

13.2 

Postcombustion  temperature  (°C) 

987 

915 

987 

897 

1085 
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sated  by  an  equivalent  loss  in  the  electrochemical  reaction, 
CO  oxidation  being  more  entropic  than  H2  oxidation  (gstack 
143  kW  versus  135  kW). 

Partial  oxidation  (POX)  reforming  of  the  biogas  obviously 
leads  to  lowered  electrical  efficiency  (43%)  due  to  partial 
fuel  consumption  already  before  entry  to  the  stack.  Cooling 
by  the  reformer,  which  approaches  autothermal  conditions, 
is  far  less  (29  kW)  and  A  has  to  be  increased  from  3  to  4.5 
to  keep  a  similar  thermal  gradient  of  approximately  200  K 
between  stack  inlet  and  outlet.  With  H2  fuel,  the  situation  is 
similar  to  the  case  of  partial  oxidation,  only  cooling  require¬ 
ments  are  yet  more  extreme  (A  =  6)  as  no  heat  is  evacu¬ 
ated  through  reforming.  For  the  hypothetical  direct  methane 
fed  system,  finally,  an  ultimate  theoretical  HHV  efficiency 
of  57%  is  found.  Heat  generation  on  the  stack  is  very  low 
(44 kW)  and  AT  maintained  at  200  K  for  an  air  excess  A 
equal  to  2. 

5.  Conclusion 

A  100  kW  class  SOFC  system  operated  on  biogas  from 
sewage  sludge  digestion  was  defined  in  a  process  flow  en¬ 
vironment.  Real  data  from  an  existing  sewage  plant  were 
used  as  input:  43m3/h  flow  of  biogas  containing  63%  CH4 
and  35%  CO2.  Using  reforming  with  added  steam  in  a  pro¬ 
portion  appropriate  to  avoid  carbon  formation  (S/C  =  0.5), 
power  generation  from  269  kW  input  methane  was  131  kW 
(48.7%  efficiency).  Fuel  conversion  was  set  to  80%  and  min¬ 
imal  cell  voltage  to  0.65  V,  requiring  2400  cells  of  20  cm  x 
20  cm  (361  cm2  active  area  each)  for  the  stack.  Cooling  of 
the  stack  was  illustrated  to  be  achievable  to  roughly  equal 
parts  by  reforming  and  cathode  air  flow  excess  (A  =  3). 
Thermal  efficiency  was  almost  40%  for  total  efficiency  of 
88.5%. 

Hot  and  cold  composite  curves  and  exergy  curves  of  the 
system  indicate  further  improvement  potential  in  internal 
heat  exchange  and  high  grade  heat  exploitation  from  stack 
and  afterburner.  This  will  be  explored  in  further  study. 

The  biogas  steam  fuel  behaved  essentially  identically  to 
steam-reformed  methane.  Using  partial  oxidation  of  biogas, 
electrical  efficiency  drops  to  under  43%  while  A  needs  to  be 
raised  to  4.5  to  maintain  a  tolerable  200  K  thermal  gradient 
over  the  stack. 
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